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Abstract 

The thermal behavior of a one-dimensional antiferromagnetic chain doped by donor impurities 
was analyzed. The ground state of such a chain corresponds to the formation of a set of fer- 
romagnetically correlated regions localized near impurities (bound magnetic polarons). At finite 
temperatures, the magnetic structure of the chain was calculated simultaneously with the wave 
function of a conduction electron bound by an impurity. The calculations were performed using an 
approximate variational method and a Monte Carlo simulation. Both these methods give similar 
results. The analysis of the temperature dependence of correlation functions for neighboring local 
spins demonstrated that the ferromagnetic correlations inside a magnetic polaron remain signifi- 
cant even above the Neel temperature T/v implying rather high stability of the magnetic polaron 
state. In the case when the electron-impurity coupling energy V is not too high (for V lower that 
the electron hopping integral t), the magnetic polaron could be depinned from impurity retaining 
its magnetic structure. Such a depinning occurs at temperatures of the order of Tjv. At even higher 
temperatures (T ~ t) magnetic polarons disappear and the chain becomes completely disordered. 



PACS numbers: 75.30.Kz, 75.30.Vn, 75.50.Pp, 64.75. +g 



I. INTRODUCTION 



The tendency to formation of inhomogeneous charge and spin states and to phase sepa- 
ration is of fundamental importance for the physics of manganites and other systems with 
strongly correlated electrons. It is now a common belief that the nature of the colossal mag- 
netoresistance effect is closely related to phase separation phenomenon The most widely 
discussed type of phase separation is the formation of small magnetic droplets (magnetic 
polarons or ferrous) as a result of the self-trapping of charge carriers in an insulating antifer- 
romagnetic matrix. This type of phase separation, first discussed by Nagaev in his analysis 
of electron states in magnetic semiconductors \2\, is now actively used for interpretation of 
experimental data on manganites at different doping levels l|, y|- 

In this connection, a problem arises, whether these self-trapped electron states can survive 
above the characteristic temperature of the magnetic ordering. The possible existence of the 
ferromagnetic correlated regions at the paramagnetic background, "temperature ferrous", 
was first discussed in detail by Krivoglaz (see also [^). Indeed, there are a lot of experi- 
mental indications (coming mostly from the susceptibility and ESR data) that manganites 
in the paramagnetic state are rather inhomogeneous. Moreover, the analysis of the high- 
temperature transport properties of manganites based on the assumption of the existence 
of ferromagnetically correlated regions in the paramagnetic phase gives a viable picture of 
the inhomogeneous state ^. 

In our paper, we study the evolution with temperature of a magnetic polaron state over 
a wide temperature range using a simple model such as an antiferromagnetic chain doped 
by donor impurities. In spite of its simplicity, this model reveals the possible existence of 
different kinds of ferrous p, Q and seems to capture the essential features of the electronic 
phase separation in manganites at low doping levels. At low doping, the impurity potential 
is a relevant parameter and we start from the state corresponding to ferrous bound by 
impurities. The density of impurities is assumed to be small, so that ferrous do not overlap, 
and the system as a whole remains insulating. We demonstrate that ferrous are rather stable 
objects and ferromagnetic correlations persist at temperatures much higher than the Neel 
temperature of the host chain. 

The structure of the paper is as follows. In Section |nj we formulate our model for the 
antiferromagnetic chain doped by donor impurities. In Section ITTH two different methods to 
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calculate the partition function for such a chain, using an approximate variational method 
and a Monte Carlo simulation, are presented. Finally, in Section HVl we discuss the charac- 
teristic features of the model in different temperature ranges. 

II. FORMULATION OF THE MODEL 

We consider a one-dimensional chain of antiferromagnetically coupled local spins. Non- 
magnetic donor impurities are distributed homogeneously with a period L (in lattice constant 
units) along the chain. The system is described by the double-exchange Hamiltonian, 

H = He\ + j'''^^Sg- Sg+l , (1) 

9 

~ ^ 4,^ 4) ^ ^, + ^imp , (2) 

Himp = -^Vgal,^ag,s- (3) 

9,s 

In Eqs. (pQ) and Sg is the spin of the magnetic ion located at site g, treated as a classical 
vector, symbols a^^, ag^g denote the creation and annihilation operators for the conduction 
electron with spin projection s at site g, and a are Pauli matrices. The second term in Eq. 
is the antiferromagnetic exchange between local spins. The two terms in i^ei describe the 
kinetic energy of conduction electrons bounded by impurities, which are located between 
sites with g = iL and g = iL + 1 (z is an integer), and the Hund's-rule coupling between 
the conduction electrons and the local spins, -ffimp is the electrostatic interaction between 
an impurity and conduction electron. Parameters V^, A, t, and J' of Hamiltonian (0) are 
considered to be positive. The hierarchy of this parameters is A > t,Vg ^ J', as usual in 
the double exchange approximation. 

As we pointed out in Introduction, the electron-impurity interaction term is needed to 
describe the low-doping limit of the model, when the material is insulating. In our previous 
discussions on this model [3,0, this potential was chosen as a deep square well. Although 
this is a good approximation for T = 0, here we prefer to use a more realistic form for this 
interaction, namely the Coulomb potential. This allows us to find self-consistently the wave 
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function, and hence the size of ferron, which becomes temperature-dependent. In fact, the 
calculations show that results are not strongly affected by the specific form of the impurity 
potential. 

We consider the low-doping limit of the model, which corresponds to large L. This allows 
us to neglect the interaction between electrons corresponding to different impurities, and 
also to restrict our calculation to a finite portion of the chain. In the part of the chain under 
study, of length L, there just one impurity located in its middle between two magnetic sites. 

Hamiltonian ()H2|) is rotationally invariant, so we can choose angles Vg between two neigh- 
boring spins in the chain as variables characterizing the magnetic structure of the local spins. 
In the limit A — oo, only one of the spin components contributes to the low- lying states. 
Starting from operators aj^, we make a transformation to operators cj, which describe the 
creation of a conduction electron with its spin directed along each local spin Sg (spinless 
fermions). This modifies the hopping term that now depends on the magnetic structure of 
local spins. Using angles Vg, we can write Hamiltonian as: 

L/2-2 

H = Hei + J COS Ug , 

g=-L/2 

L/2-2 L/2-1 t 

^ei = -t (4+1^9 + i^9+i) cos^-V J2 \a + i/2\ ' 

9=-L/2 g=-L/2 1^ / I 

where the Coulomb form of the impurity potential is introduced explicitly, and J = J'S"^. 

III. CALCULATION OF THE PARTITION FUNCTION 

Using Hamiltonian (^, we can write the partition function in the following form: 

Z = I / sin UgC-f^J ^9 ] I e-/^^ei ) I (5) 

n=0 [g=-L/2 \{ J J 

where f3 = 1/T, and {ipn) are the basis electron wave functions, which are assumed to be 
orthogonal: 

L/2-1 L/2-1 

1^-) = E ' E ^fVf =5n,^- (6) 

g=-L/2 g=-L/2 
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In this expression, we suppose that il)g are real numbers, because the matrix elements of 
operator ifgi in the |0) basis are real. Note that we can use any appropriate basis to 
calculate the partition function. If Itpn) are the eigenfunctions of the operator ifeb the last 
factor in Eq. (jSj) can be written in the form: 

I e-^^-i ) = e-^ I 1^" ) = e-P^- . (7) 

In this case, wave functions Itpn ) depend, of course, on the canting angles. 



A. Approximate partition function 

Partition function © can not be calculated analytically. In this subsection, we propose an 
approximate approach to the evaluation of the partition function and physical observables. 
The main point of this approach is the substitution in Eq. (0): 

{^n I ) ^ I 1^" ) , (8) 

where effective wave functions | do not depend on canting angles Vg. The procedure of 
finding these wave functions is described below. The accuracy of such a procedure is tested 
by the Monte Carlo simulations. The physical meaning of this substitution is that the main 
contribution to Z comes from the term corresponding to the ground state of the conduction 
electron. Since each electron in the ground state is trapped by its impurity, we can conclude 
that the its wave function only slightly depends on the canting angles. 

After substitution of Eq. (jHl) into Eq. the approximate partition function can be 
written in the form: 



L-l 



^appr _ ^exp 



n=0 



g=-L/2 



in) 



9 + ^21 



X 




duq sin z/g exp 



"/?(2t?l:\V^i"^os^-Jcos.^ 




(9) 



where we use Eq. (jD) for operator H^i- 

— (n) 

The wave functions ifjg are found by minimization of the approximate free energy 
^appr _ _rp ^■appr -^ith rcspcct to ip g . Howcvcr, we should take into account that wave 



functions ^'g are not independent because of the orthogonahty conditions (0). Calculating 
the derivatives dF'^^^'^ / difj^ , we obtain the following system of nonlinear equations: 

n 

((^el)(n,„)V^^"^) - V Krn^t^ = , n = 0, 1, . . . L - 1 , (10) 

^ ^9 — n 

m=\j 

where H^x is the matrix of the operator i^ei in the cj|0) basis , A„ m (m < n) are the 
Lagrange multipliers, and symbol (. . . )(n,m) denotes the following averaging procedure: 
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dugsmugf{ug)e 



tyi^g+li^g +i^g+li^g j COSy - JCOSZ/^ 
{fi^g))(n,m) ^ „ r /— fnl — Cml —(m.) —(n)\ Vr. I 1 ' ^^^^ 



dUg sin UgC 



P t (^g^l^g""^ + ^1+1^1"^) COS ^ - J COS Z/g 







The system of equations (fTn|) with additional conditions (jHl) is solved numerically. Note that 
the approximate free energy satisfies the inequality F^pp"^ > F, because for any canting angle 
the following condition is met 

(^J i/el|?J >i?n . (12) 

Let us make some remarks concerning Eqs. ()10|1 and (fTT|) . It can be easily seen from 
these equations that parameter A„ is the mean energy of the electron in the nth state, 
that is An^n = (V'n I (Hei) {n,n) \^n) ■ Wc cau Say that \tpn) are the states of an electron 
moving in the averaged local spin background corresponding to these states. But, strictly 
speaking, only the state with n = ("ground" state) is the eigenstate of the averaged 
operator (ifei)(o,o)- 13 clear that in the case when the difference between energies A„_„ 
is large enough in comparison with temperature, that is /3(A„4.i „+i — A„ „) ^ 1, we can 
neglect contributions to the free energy from all states with n > 0. 

The approximate approach described here allows us to calculate also the physical observ- 
ables. Let us consider first the local spin sector of the problem. The mean value of the 
function F{{9g, cpg}) can be easily calculated only in the special case when: 1) F depends 
only on Ug, and 2) it is possible to separate the variables in the {L — l)-dimensional integrals 
over the canting angles. The simplest case is: 

L/2-2 

F{{eg,<Pg})= f^^S^- (13) 

g=-L/2 



It is useful to introduce the notation: 

/ L/2-1 -An)—{rn) 

Znm = exp (3V V ^ , I X 

dvg sin z/g exp 
The mean value of fivg) then reads: 




(^t[^g+1^9 +^3+1^9 j COS y - /3 J COS Z/^ 



.(14) 



1 



(15) 



n=0 

To calculate the mean value of the operator O corresponding to the physical quantity in 
the electron sector, we should made the transformation: 

L-1 

(16) 



m,=0 



where On,m = yPn I ^ |V'm) the matrix elements of the operator O in the basis 
The mean value of O can be written in the following form: 
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^■appr 



(17) 



n,m=0 



We can also find mean values of physical observables in more complicated cases. For 
example, the approximate value of the mean energy of the system is: 



^■appr 



1 

ippr 



■'m,n 



n,m=0 



L/2-2 — (n)— (m) 



g=-L/2 



L/2-2 



l^7 + l/2| 



COS 



g=-L/2 

J 

^■appr 



2 / (n,m) 



L-1 L/2-2 



J2 ^n,"(cOSi/3)(„_„) . 



(18) 



n=0 g=-L/2 

At relatively low temperatures, /3(A„+i^„+i — A^^^) ^ 1, formulas (fTHjl . (fTTjl . and (fTR|l for 
mean values have a simpler form. For example, the mean value of f{i^g) is equal to (/(t'g))(-Q g)- 
To test the approximate approach described above, we calculate also the partition function 
and physical quantities characterizing the temperature behavior of the bound magnetic 
polaron using Monte Carlo simulations. 



B. Monte Carlo simulations 

Partition function (0) can be numerically calculated using standard diagonalization sub- 
routines and a Monte Carlo (MC) simulation. For a given local spin configuration, the 
electronic Hamiltonian is quadratic in the Cg operators and can be exactly diagonalized. 
The remaining integrals over the local spin configurations can be calculated using a classical 
MC simulation. The MC simulation of the local spin system is done using the Metropolis 
algorithm (see, for example, Ref. j^). The Metropolis algorithm generates a set of local spin 
configurations according to some probability distribution for the local spin angles. We use: 

L/2-2 L-1 

P{{ug})= H (^e-^'^cos'^5)^e-^^«({^9}) . (19) 

g=-L/2 n=0 

Additional weight factor coming from the volume elements of integrals (0) is included in 
the course of calculating the physical observables. Note that for the diagonalization of the 
electronic Hamiltonian, one needs to know the local spin configuration, and for the MC 
simulation of the local spin system, one needs to know the ground state of the electronic 
Hamiltonian. We proceed as follows. Starting from a given spin configuration, we diagonalize 
the electronic Hamiltonian, and calculate P {{i^g})- Then we choose a site g at random, and 
temporarily change Ug u'g. The electronic Hamiltonian is again diagonalized and the 
new P ({t'g}) calculated. The change in canting angles is finally accepted according to the 
Metropolis algorithm. We define a MC step as a number of these spin reorientations equal 
to the number of canting angles, L — 1. Typically, we use series of several thousands of MC 
steps for each set of parameters of the Hamiltonian and T. A half of these steps is used 
to thermalize the system, and another half is used to calculate the physical observables. 
Between each two MC steps, we discard 10 extra MC steps to avoid correlations between 
consecutive measurements. For the spin reorientations, the new value of the corresponding 
canting angle is randomly chosen. Open boundary conditions are used at the ends of the 
chain. Note that in contrast to the approximate method described in the previous subsection, 
the MC simulation can reveal some effects related to the dynamics of the conduction electron, 
the time being defined by the MC steps. 
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IV. RESULTS 



We calculate partition function of the model (0} for the chain with 16 sites at different 
values of parameters V, t and J . This length is representative of the low doping limit, and 
it is found that further increase of it does not change the obtained results. The impurity 
is placed in the middle of the chain between sites with g = —1 and g = 0. To compare 
results given by variational method and MC simulation, we calculated the mean energy of 
the system E in the temperature range < T < 2.5J. The calculations of function E{T) 
for different parameters of the model show that difference between both methods does not 
exceed several percent. Some test MC calculations at higher temperatures also demonstrate 
a good agrement with approximate procedure. So, we believe that our variational approach 
is a good approximation even at high temperatures, and calculate physical quantities at 
T > 2.5 J using this method only. 

As was already mentioned in Section m our methods allow us to calculate the wave 
function for the conduction electrons. From this, we can estimate the size of the bound 
ferron. We have used different set of parameters of the Hamiltonian, and found that the 
size of magnetic polaron depends only slightly on the interaction between impurity and 
conduction electron, V. For example, for V = lOOJ magnetic polaron contains 4 magnetic 
sites, whereas for V = lOJ it extends over 6 sites. 

To analyze the evolution of the magnetic polaron state with the growth of temperature, 
we calculate the correlation functions of neighbouring spins in the chain, i.e. mean values 
{Sg- Sg+i) = S^{cosh'g). FoT thc temperature range < T < 2.5 J, we have performed the 
calculation of correlations function using both methods described above. Corresponding 
curves are shown in Fig. ^ for V = 100 J, t = 50 J. It can be seen that both methods 
give similar results. Figure ^ shows that for the four sites closest to the impurity, we have 
ferromagnetic correlations, i.e. (cosz/g) > for g = —2,-1,0. The remaining part of the 
chain exhibits antiferromagnetic correlations. At zero temperature, local spins inside the 
ferron are nearly parallel each other, i.e. (cosug) ^ 1, whereas neighboring spins in the 
remaining part of the chain are antiparallel, i.e. (cosz/^) ^ —1. 

Due to one- dimensionality of the model, the magnetic phase transitions in the chain do 
not occur. Nevertheless, we can introduce some characteristic temperatures describing our 
system. The antiferromagnetic correlation functions, (cosug), between neighboring spins far 
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FIG. 1: Correlation functions of local neighboring spins in a finite 16-site chain, {Sg- Sg^i), versus 
temperature in the range < T < 2.5J. The values of local spins are normalized to 5" = 1. The 
parameters of the Hamiltonian are V = 100 J, t = 50 J. The panels corresponds to the different 
sites (a) g = —1, (b) g = 0, (c) g = 1, and (d) g = 6- The plots for sites g = 2, . . . ,6 are nearly the 
same. Also notice that the chain is symmetric with respect to the impurity located at g = —1/2. 
In all panels, solid lines correspond to the variational method described in the text, and points to 
the MC simulations. 

from the impurity exhibit similar behavior. We may define temperature T)^p corresponding 
to the Neel temperature of 3D Heisenberg antiferromagnet in a conventional way, as the 
point of the steepest change in this correlation function (that is the point where its curvature 
has a maximum). This temperature does not depend on V, and equals approximately to 
T^F = 0.56 J, which is close to the mean-field estimate of this parameter. From Fig. ^ we 
see that it corresponds to the point where (cos z/g) |t=t*j^ ~ —1/2. 

The plots in Fig. ^ demonstrate that the ferromagnetic correlations are significant for 
the first and second magnetic neighbors to the impurity, even at high temperature. The 
AF correlations in the rest part of the chain decay much faster with temperatures. So, the 
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ferron is a stable object that does not disappear even at T T\p. 

Let study the temperature range in which the ferromagnetic correlations start to decay 
steeply (see Fig.|2I). We introduce, similarly to T^^, a set of temperatures characterizing the 
magnetic polaron: T^*, Tg*, . . . , T^, n being the last site showing ferromagnetic correlations at 
zero temperature. Namely, is the temperature corresponding to the maximum curvature 
of the plot shown in Fig. |21for the correlation function for the sites nearest to the impurity. 
From Fig. 121 we see that this temperature corresponds to the value (cos z/_i) |t=t* ~ +1/2. 
Similarly, Tg* corresponds to the value at which (cos z/q) |t=T2* ~ +1/2, and so on. In the 
case of the sets of parameters corresponding to Fig. |2l n = 2 for V = 100 J, and n = 3 for 
V = 5J. These temperatures are T; = 1.45J, T* = 8.34J for V = lOOJ, t = 50J, and 
T* = J, T* = 3.08J, T* = 3.92J for V = 5J,t = 50J. At temperature T = in the case 
of V = lOOJ, and at T = Tg* for V = 5J, the ferromagnetic correlations inside the ferron 
start to decay, and at T > T^, the ferron state completely disappears. 

The calculation shows that the ferron state is more stable at larger values of V. How- 
ever, this dependence is not very much pronounced: the change in the characteristic tem- 
peratures is by a factor of 2, when V changes from lOOJ to 5J. The temperature Tj*, 
which characterizes the complete decay of the ferron, is in the range of 0.1 — 0.2t. The 
value of can be obtained on the basis of following simple estimates. If we note that 
at Tj* all correlation functions except the first one go to zero, and that the wave func- 
tion of the electron does not change very much with temperature, we can approximate the 
gain in energy of having a ferron state as (^ t-i/'^V's+i cos (z/g/2)). Using that at T = T^, 
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(cosz^_i) ~ 1/2, and assuming that the wave function is uniformly distributed over the size 
of ferron, we have \ipg\'^ = 1/4 for the 4-site ferron, {ipgl"^ = 1/6 for the 6-site ferron. Then, 
~ (^l^oP cos (z/_i/2)) ^ 0.1 1 — 0.2 1. Also note that the relationship between Tj* and t is 
similar to the relationship between the Curie temperature Tq and t for the double exchange 
model. This similarity could have the same physical nature since in the latter model, the 



ferromagnetic ordering is related to t 
conduction electrons on a lattice site 



le same hoping integral multiplied by the density of 
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Using our approximate variational method, we can analyze the case of —>■ 0. The 
plots for the correlation functions are very similar to those in of Fig. ^p. In this limit, we 
have n = 4, so the ferron contains 8 sites. The values of the corresponding characteristic 
temperatures are T* = 2.97J, T* = 2.57J, T* = 1.31J, and T* = 0.14J. Therefore, in our 
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FIG. 2: Correlation functions {Sg- Sg+i) = (cosfg) versus temperature at parameters V = 100 J, 
t = 50 J (a), and V = 5 J, t = 50 J (b). In both cases, the upper curve corresponds to the correlation 
between spins at sites near to the impurity (cosi^_i), the next one is the correlation between next 
two spins (cos uq) and so on. The lower curve is the correlation between last two spins in the chain 
(cosine)- In panel (a), ferromagnetic correlations exist for first two sites near to the impurity, and 
magnetic polaron contains 4 sites, whereas at smaller V (panel (b)), it contains 6 sites. 

approach, we can describe both the ferron bound by impurity, and depinned from it. Some 
indication of this possible depinning at small values of V could come from the presence of 
small ferromagnetic correlation corresponding to the fourth pair at high temperatures (see 
Fig.Eb). 

To study the possible depinning in a more detailed way, we calculate the quantum- 
mechanical mean position g = Z~^^ J2J2{'^9^^\ 9 "^g^^) of conduction electron for 

n g 
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FIG. 3: Probability distribution of the quantum-mechanical mean position of the conduction 
electron, g, for different values of temperature, and V = 5J. The probability distribution broadens 

with the growth of temperature. The insets show the value of ^ as a function of the MC time. 
At T > 1.5 J there is a large scatter of g, which is the signature of the ferron depinning. Vertical 
dot-dashed line indicates the impurity position. 
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FIG. 4: Standard deviation of g, Ag, as a function of T. Different curves correspond to different 
impurity potentials: V = 0.1, 1, 5, and 100 J (corresponding the first to the highest, and the last 
to the lowest curve, respectively). 

each MC time, where Z^i = Dependence of mean position g on MC time is 

n 

shown in insets to Fig. 01 In this figure, we also show the histograms illustrating probability 
distribution of g in the temperature range < T < 2.5J, and for the set of parameters 
V = 5J,t = 50 J. In the insets we can see that at temperatures T > 1.5 J, there appears a 
large scatter in the values of mean position g. Simultaneously, the probability distribution 
broadens. This means that an electron becomes depinned and can move far away from 
the impurity. Since the depinning temperature is, at least, lower than the temperatures 
T^, T2 for the same values of V, and t, we can conclude that the ferron retains its magnetic 
structure after depinning. At sufficiently large values of V, the depinning of a ferron is not 
observed in the temperature range under study. 

In Fig. we plot the standard deviation of the histograms of (7 as a function of temper- 
ature for several values of V. In all cases, the standard deviation has a maximum for some 
7"max > J- At relatively low temperatures, Ag grows due to the increase in the thermal 
fluctuations of the ferron position. At higher temperatures, the ferromagnetic correlations 
forming the ferron gradually decay, and the electron undergoes a kind of Anderson localiza- 
tion in disordered spin background (the situation typical of the double exchange mechanism). 
Therefore, Ag eventually decreases independently of the value of Coulomb potential V. 

We can also find the depinning line using the variational method. Let us consider the 
ferron state bound to the impurity and the state at 0. The latter case corresponds 

also to the free ferron located far from the impurity. The energy difference between bound 



14 



32 












16 
8 




I 




dep 


II 


4 
1 






rn* 
^1 




' ~ " - ^ 2 

n = 2 








- - ^ - ^ ^7^3* 


0.5 












0.25 












0.125 
0.0625 




'/.^ Tl 
= 8 ~ - - 




n = 6 \ 


n = 4 >^ 



0.0625 0.25 1 4 16 64 256 

VIJ 



FIG. 5: Phase diagram in the V — T plane for a chain with a low density of bound ferrons. The 
value t = 50 J is fixed. Solid line corresponds to T*, that is, to the complete decay of the ferron. 
Dashed lines are the T*!^ plots, depicting the decay of ferromagnetic correlations at the end pair 
of spins in the ferron. Horizontal line corresponds to the temperature T'^p. Dot-dashed line is 
the depinning line obtained by estimations described in the text. Above this line, ferrons are not 
bound by the impurities. Regions I and II above the "melting" line T* correspond to the free 
electron in the chain, and to the electron bound by the impurity in the absence of the ferron state. 

and free ferron states then equals to E{0,Tdp) — E(y,T). The depinning temperature T^p 
can be estimated as: 

E{0, Tap) - E{V, Tdp) = T^p , (20) 

where the mean energy of the system E{V,T) is calculated by variational method using 
formula (fTHj) . Note that the depinning temperature calculated in such a way turns out to 
be in several times greater than its value estimated by Monte Carlo simulations, because 
in latter case we consider the states corresponding to the ferron located not far from the 
impurity. 

In Fig. we present the phase diagram of the chain in — T plane. The tempera- 
tures Tj*,T2*, ■ ■ ■ ,T* as a functions of the impurity potential are shown. We also plot the 
temperature T^^p and the depinning line calculated by the variational method. 
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V. CONCLUSIONS 



We show that the simple model of an antiferromagnetic chain with impurities captures 
the essential features of the phase separation in low-doped manganites. We find a set 
of temperatures characterizing the thermal evolution of the magnetic structure of such a 
chain. We demonstrate that the magnetic state of the chain is rather inhomogeneous and it 
can be described by a set of magnetic polarons. These magnetic polarons are rather stable 
objects existing at temperature even much higher than the Neel temperature of the undoped 
chain. With the growth of temperature, ferromagnetic correlations inside the ferron start 
to disappear in a step-wise way, up to some temperature T* ~ O.lt — 0.2 1, and then ferron 
state is no longer stable. The dynamics of ferrous depends on the magnitude of the electron- 
impurity coupling, V. At high values of V, ferron remains tightly bound by impurity, while 
at lower V, the ferron can be depinned from impurity at a characteristic temperature lower 
than T^*, and therefore retains its magnetic structure after depinning. These results can 
be summarized in the phase diagram in the V — T plane presented in Fig. El We can see 
that in region below the rather flat curve versus V, there are the regions of stability of 
different kinds of ferrous, from n = 2 to 8 sites. In each such region, we show the dashed 
line corresponding to the melting of the last ferromagnetically correlated pair in the ferron. 
So, between the solid and each dashed line, we have the region where the ferron gradually 
decays. We also show the depinning line calculated by variational method as it was discussed 
at the end of the previous section. We see that there is a large region between the depinning 
line and T* where ferron is depinned. We have extended this depinning line above the Tj* 
curve to show the possible domains of existence of bare electron bound to the impurity and 
freely moving along the chain. 
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